Non-Hermitian Luttinger Liquids and Vortex Physics 
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As a model of two thermally excited flux liquids connected by a weak link, we study the effect of a 
single line defect on vortex filaments oriented parallel to the surface of a thin planar superconductor. 
When the applied field is tilted relative to the line defect, the physics is described by a nonhermitian 
Luttinger liquid of interacting quantum bosons in one spatial dimension with a point defect. We 
analyze this problem using a combination of analytic and numerical density matrix renormalization 
group methods, uncovering a delicate interplay between enhancement of pinning due to Luttinger 
liquid effects and depinning due to the tilted magnetic field. Interactions dramatically improve the 
ability of a single columnar pin to suppress vortex tilt when the Luttinger liquid parameter g < 1. 
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INTRODUCTION 

The past decade has seen considerable work on the 
statistical mechanics and dynamics of thermally excited 
vortices in type II high-temperature superconductors . 
The competition between interactions, pinning and ther- 
mal fluctuations gives rise to a wide range of novel phe- 
nomena, including a first order melting transition of the 
Abrikosov flux lattice into an entangled liquid of vor- 
tex filaments and a low-temperature Bose glass phase 
with vortices strongly pinned to columnar defects Q ■ 

A convenient way of understanding interacting flux 
lines is provided by the formal mapping between the 
classical statistical mechanics of (d+l)-dimensional di- 
rected flux lines and the nonrelativistic quantum mechan- 
ics of d-dimensional bosons. In this mapping, flux lines 
traversing the sample along the direction of the external 
magnetic field H = Hz correspond to boson world lines 
propagating in imaginary time r. The classical partition 
function of thermally excited vortex lines is proportional 
to a quantum mechanical matrix element. The thickness 
of the sample in the ^-direction, L z , corresponds to the 
inverse temperature (3% of the bosons, while thermal fluc- 
tuations of the vortices, due to finite HbT, play the role 
of quantum fluctuations of the bosons, controlled by h~. 

If the direction of the external magnetic field does not 
coincide with z, the direction of the columnar defects, 
it is convenient to separate the transverse component of 
the field H± from the parallel one H\\ along z. When 
H± <C -ffy, the transverse component H± plays the role 
of a constant imaginary vector potential for the bosons 
0,3- The corresponding quantum Hamiltonian is non- 
hermitian, with new and interesting properties. Stimu- 
lated by vortex physics, there has been considerable in- 
terest in models of noninteracting quantum particles in 
a constant imaginary vector potential with a disordered 
pinning potential. 

Less is known about non-hermitian models with inter- 
actions. A disordered array of parallel columnar defects 
leads to a strongly pinned low-temperature Bose glass 



phase. For H± less than a critical value H±, this phase 
exhibits a "transverse Meissner effect" , such that the vor- 
tex filaments remain pinned parallel to the columns even 
though the external field is tilted away from the column 
direction. 

In this paper we study the effect of a single columnar 
pin (or an equivalent linear defect) on the statistical me- 
chanics of thermally fluctuating vortex lines confined in a 
thin, superconducting slab (see Fig. QJ. With exception 
of a brief discussion in Q and the analysis of Ref. 
little has been done on vortex physics in the limit of a 
dilute concentration of columns (or twin planes). If the 
average spacing between defects is d, this is the regime 
Hh ^> £?$ = $o/^ 2 where £?$ is the "matching field" and 
$0 is the flux quantum. The large number of "intersti- 
tial" vortices between defects could be locally crystalline 
or melted into a flux liquid. A third possibility (perhaps 
the situation most closely related to our calculations in 
1+1 dimensions) is vortices in a "supersolid" phase, with 
extended correlations in both the translational and boson 
order parameters 0. 
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FIG. 1: Schematic snapshot of vortex lines and a columnar 
pin. The near-vertical segments occur at locations of maxima 
of the average density. For finite tilt these maxima become 
negligible for \x\ > fx. Here, the z-axis is denoted by t. 

The feasibility of studying vortex physics in sam- 
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pies which are effectively (l+l)-dimensional was demon- 
strated by Bolle et al. [3 in thin samples of NbSe2, where 
the effect of point disorder on interacting vortices near 
H c i was observed. Similar experiments might be possible 
on thin high-T c samples where columnar defects could be 
implemented mechanically by cutting a thin "notch" . A 
related problem in (2+1) dimensions concerns the effect 
of a single twin plane or grain boundary on vortex mat- 
ter, where point disorder may lead to algebraic decay of 
density correlations of the Abrikosov flux lattice sim- 
ilar to those discussed below. A single such plane has a 
similar pinning effect on bulk flux lines as does a linear 
defect in (l+l)-dimensional systems. 

In the following we will consider a single columnar de- 
fect in a system of interacting flux lines in (1+1) dimen- 
sions. Although results can also be obtained using a flux- 
line related phonon formalism [fij , we here found it conve- 
nient to work with an equivalent quantum Hamiltonian 

us 



Hamiltonian Q 



H = - 
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where V(|a;|) is a short-range repulsive vortex interaction 
potential, ^(x) annihilates a bosonic flux line, n(x) = 
^(x) 1 i(x) is the boson number density and m is the 
vortex tilt modulus. The imaginary vector potential h = 
&qH±/ (4ttT) arises due to the tilted magnetic field and 
to is the strength of the defect modelled by a J-potential 
at the origin. In the following we set ksT — 1 (i.e. K = 1 
in the quantum model). 

Without the local potential and the non-hermitian 
term this model has been well studied 0, . In par- 
ticular, Haldane has shown that this spinless Lut- 
tinger liquid exhibits a line of critical points with con- 
tinuously varying exponents. His calculation is based on 
the bosonization technique, where the boson field 
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dx 
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^i2-Krn(riQX-\-u(x)) i(j>{x) 



(2) 



is represented in terms of a boson phase operator 4>{x) 
and (dimensionless) phonon operator u{x). The two 
fields satisfy the commutation relation [<fi(x),u(y)] = 
(i/2) sgn(cc — y). With applications to vortex physics 
in mind, we have extended the bosonization approach 
and work on quantum impurities |l2l to the non- 
hermitian case h > and calculated asymptotic low- 
energy properties for the model lfl|. 

In addition, we have performed a non-perturbative nu- 
merical analysis using the Density-Matrix Renormaliza- 
tion Group (DMRG) for a discretized version of the 
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corresponding to a nonhermitean Bose-Hubbard model 
where rn = b\b i and the hopping is t = l/2m (for unit 
lattice constant). In the following we set t = 1. We work 
in the canonical ensemble, fixing the density of bosons per 
site no- We have retained an onsite and a next-neighbor 
interaction, which turn out to be sufficient to qualita- 
tively describe the full phase diagram. Furthermore, for 
computational purposes we allow at most 2 bosons per 
site, which effectively renormalizes the on-site repulsion. 
The lattice model is a good approximation to 
for small filling uq (average number of bosons per site). 
Our calculation is based on an extension of the DMRG 
to non-hermitian systems with complex eigenvalues and 
eigenvectors (for details see 0). 

In the hermitian case h = without impurity we have 
first calculated the Luttinger liquid parameter g which 
governs the long-wavelength behavior of correlation func- 
tions and is important to understand the response at 
finite tilt. We adapt the DMRG work of Ref. [IcJ to 
periodic boundary conditions, essential for the study of 
persistent currents (i.e. arrays of tilted vortex lines) dis- 
cussed below. We focus exclusively on the superfluid 
(Luttinger liquid) phase. Via DMRG we have calculated 
the boson correlation function, which from conformal 
field theory is expected to behave as < ^(x)^(0) >^ 
|Lsin(7ra;/L)|~ 1 / 29 . We have verified this behavior nu- 
merically with high accuracy and have extracted g by a 
fit to the data (see Fig. EJ). We also determined g from 
the compressibility and the finite-size dependence of the 
ground state energy |17| with excellent agreement be- 
tween the values of g obtained by both methods. For 
arbitrary short range potentials (in continuum or lattice 
models) we have derived the general low-density result 
jkI - 2ano + 0{a 2 nfy where a is the two-particle scat- 
tering length. For our lattice Hamiltonian J3J we find 
o= -(8t 2 ~4tV-UV)/(2tU + UV + AtV). As shown in 
Fig. [2 this asymptotic result is in good agreement with 
the numerical data. 

We now include the pinning term proportional to eo 
in Eq. (Q. In order to determine the relevance of this 
term at long wavelengths we have performed a perturba- 
tive renormalization group (RG) analysis. We obtain the 
following renormalization flow of the pinning strength 



co(0 = eo(Jo) ( j 
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(4) 



where I is an effective length scale or inverse cutoff mo- 
mentum. For g > 1 the renormalized coupling to flows 
to zero at long length scales while for g < 1 it diverges. 
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FIG. 2: Luttinger-liquid parameter g vs. density. The 
dashed lines show the analytic result at low densities. 




FIG. 3: Friedel oscillations of the flux line density vs. dis- 
tance x from the defect, calculated numerically for U = 10, 
V — 4, L — 128 and no = 0.25, corresponding to g ~ 0.72. 



Remarkably, while in fermionic systems with (generic) 
repulsive interactions one always has g < 1, the bosonic 
Luttinger liquid studied here can be tuned to either 
regime. This can be easily seen by setting U = oo, V = 
in (5J : Since hardcore bosons in Id are equivalent to non- 
interacting spinless fermions, we obtain g — 1. Smaller 
U increases g from 1 while additional next-neighbor in- 
teractions V > lead to g < 1. In the following we will 
denote the special situation g = 1 as the free fermion 
limit, for which we have replaced the DMRG by compu- 
tationally less expensive exact diagonalization. 

The irrrelevance/relevance of the pin can be clearly 
observed in the Friedel oscillations of the boson density 
Afi(i) = (n(x)) —riQ, for which we find the analytic result 



(5) 



with an interaction-dependent exponent a = | 2 g-i } f° r 

and a decay length £j_ ~ l/h. At zero tilt h we 

therefore find a pure power-law decay of the density per- 
turbation which is more rapid when the pin is irrelevant. 
When h > 0, the Friedel oscillations acquire an additional 
exponential decay. 

As illustrated in Fig. [21 the density oscillates with a 
phase set by the impurity position, and an algebraic en- 
velope before exponential decay sets in for x > In 
the vortex picture (see Fig.^), configurations are domi- 
nated by parallel, tilted flux lines at distances larger than 
£±{h) from the pin. Closer in, vortices attempt to align 
with the maxima in the density oscillations present when 
h = 0. This alignment is limited by interactions as vor- 
tices enter and leave the aligned region with increasing 
imaginary time r. The resulting vortex configurations re- 
semble a symmetric traffic jam, with vortices queuing up 
(and occasionally changing places) in the vicinity of the 
columnar defect. With our conventions, the slope of the 
lines far from the pin is h/m, so new vortices enter the 
jam at imaginary time intervals £|| ~ m/hriQ, where no is 



the linear density of "bosons" . If c is the Luttinger liquid 
velocity we expect that £_l(/i) ~ c£|| ~ 1/gh, a diverging 
length scale we confirm with our analytic calculations. 
Note that the pinning strength is reduced dramatically 
for length scales x > £,±(h) even for g < 1. 

When h > 0, the non-hermiticity leads to a finite per- 
sistent current J}, — —j^d(H) /dh in the ground state. 
This current is purely imaginary and corresponds to the 
transverse magnetization in the original flux line system 



Mi - $oImJ b . 



(6) 



The defect reduces this current, due to flux lines pinned 
even in the presence of a tilted magnetic field. Al- 
though a single pin cannot modify the bulk current in 
the thermodynamic limit £ — » oo, it creates nontrivial 
finite-size effects. Since ImJb = hN^/mL in the ab- 
sence of pinning (where Nf, = hqL), it is convenient 
to define a "pinning number" N p for vortices given by 
ImJf, = h(Nb — N p )/mL. Because Re J;, = 0, this can be 
written 



N p = N b [J b (0) - J b {e o )}/J b (0) 



(7) 



The quantity N p may be readily calculated for the free 
fermion case g = 1 where the groundstate energy is de- 
termined by filling up all the states below the "Fermi 
surface" . The pinning number obtained in this way has 
the asymptotic behaviour (L — » oo) 



N„ 



(me << h) 



(me Q ) 2 

-^ln(|e |/no), (me »h), 
h 



(8) 



results valid provided h <C ?io = N b / L. Remarkably, N p 
diverges as h — > 0. When the pinning is strong, the func- 
tional form (JHJ can be understood in terms of the aligned 
local density wave which extends out to a distance 
The Np « £,±riQ « no/h vortices entrained in this "traf- 
fic jam" do not contribute to the current. To check this 
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FIG. 4: Pinning number in the free fermion limit g = 1 for 
L = 100 and eo = 2. Note the "step" due to single-vortex 
depinning and the strong enhacement at small tilt h. 

divergence, we have also calculated the pinning number 
numerically within the lattice model Q, both in the free 
fermion limit (U = oo, V = 0) and for general inter- 
actions. Results are shown in Fig. A clearly visible 
feature is the "step" at intermediate tilt for low boson 
densities, correspondingto the single-vortex depinning 
transition at h c w meo |4j. Most prominent, however, is 
the dramatic increase in the number of pinned vortices 
at small tilt. We find similar results with DMRG for 
g < 1. In the linear response limit (hL — > 0) we find 
more generally that N p ~ L 3 ~ 2g for an irrelevant defect 
(g > 1), while in the relevant case almost all vortices 
are pinned, i.e. N p — > iV& for large system size L. The 
residual current for a relevant pin has the linear response 
form Jf,(h)\h->o ~ hL 1-1 / 9 , which vanishes as L — > oo. 
A similar result for fermions was obtained by Gogolin 
and Prokov'ev [18] in the case of a real vector potential. 
We have verified this finite-size dependence of the cur- 
rent with good accuracy in the DMRG calculation (see 
Fig.EJ. 

The equivalence to a real vector potential breaks down 
at finite h. Our results are consistent with two qual- 
itatively different types of behaviour depending on the 
value of g. While for an irrelevant pin (g > 1) a simple 
power-law scaling ansatz N p (h) = h~ 3+2g §(hL) works, 
our analytic work suggests a nontrivial logarithmic cor- 
rection for g < 1: 

N p (h) = (n /h) (- - 1) ln(h) + const.) (9) 

This equation is valid for hL 3> 1 and h <C no- The 
DMRG data (see inset of Fig. |SJ are consistent with this 
conjecture. 

In conclusion, we have studied the effect of a single 
columnar defect on a sea of interacting vortices in 1+1 
dimensions, in the presence of a tilted magnetic field. 
The physics is described in terms of the groundstate of a 
nonhermitian Luttinger liquid. Our calculations demon- 
strate that repulsive interactions can lead to a dramatic 
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FIG. 5: Main plot: Finite-size scaling of the current (DMRG 
results) for filling no = 0.25 and a relevant pin (g ~ 0.72). 
Notice the data collapse in the linear-response regime hL — > 
0. Inset: pinning number N p as a function of h for the same 
parameters. The dashed line gives the logarithmic behavior 
in Eq. @ with an offset of const. = 0.5. 



enhancement in the number of pinned flux lines for g < 1 
and thus to a strong transverse Meissner effect controlled 
by £±(h). Details of our analytic and numerical work will 
appear later [lflj . 
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